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Abstract 

We show that the spectral flow due to a moving vortex is identical to the phase 
slippage process, and conclude that its evaluation confirms the results by the Berry's 
phase calculation of the Magnus force. 
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Consider a rectilinear vortex moving in a superconductor with a small velocity v^, relative 
to the background of the crystal lattice, both the Berry's phase|jl[] and the direct total 
transverse forcep| calculations lead to the Magnus force: 

Fa/ = -L qi,h^ v^, x z . (1) 

Here ps is the superfiuid electron number density, = ±1 is the vorticity, and the vortex 
line is along the z-direction with the length L. The crystal lattice background of the super- 
conductor is held at rest in the laboratory frame to avoid further complications. Looking 
from the point of view of the electron fluid, a Fermi liquid, moving vortices cause phase 
slippages, generating a potential difference. Therefore the electron fluid also feels a force, 
corresponding to the measured electric field. This process has been analyzed by Josephson[Q 
and Anderson]^]. Recently, it has been reanalyzed by the calculation of the spectral or mo- 
mentum flow in the electron fluid 0, ^. Since this spectral flow has a sign opposite to that of 
the Magnus force, one might be tempted to conclude that it cancels the Magnus force. The 
purpose of the present letter is to show that those two views of looking at the consequences 
of a moving vortex are equivalent as the action-reaction forces. This equivalence must have 
been implicitly contained in the literature, but it has not been explicitly stated yet. Instead, 
some confusion still exists. [0, || 

We start by studying the momentum flow in the electron fluid due to a moving vortex 
trapped by a potential located at tq. The trapping potential can be arbitrarily weak, and 
just defines the vortex position. The electron fluid system is homogeneous other than this 
trapping potential, and the system of vortex-electron fluid is translational invariant. The 
momentum flow in the electron fluid due to the moving of the trapping potential, the moving 
of the vortex, is 

jP{t) = TrimP] . (2) 
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Here the density matrix for the whole electron fluid is governed by the equation ihp{t) — 
[H,p{t)], with P the total momentum operator and H the truncated Hamiltonian of the 
system: the lattice phonon induced effective electron-electron interaction is treated as an 
attractive one, which gives arise to the superconductivity. One may include the lattice 
dynamics into the formulation. In this case, because the electron-lattice system consists of 
a clean superconductor, our following conclusion will remain unmodified. We first evaluate 
the momentum flow from a 'safe' distance far away from the vortex core. Later we point out 
that it can be evaluated near the vortex core by the spectral flow method and the results of 
the two calculations are equal as guaranteed by the momentum conservation law. We note 
that for a slow moving vortex, the adiabatic condition holds: At a given time the density 
matrix p of the electron fluid system can be approximated by its instantaneous equilibrium 
density matrix po and the deviation p^. p = po + pi- Here pi is determined by the equation 
[H, pi] — ihpo, and 

Po(t) = E/nl^n >< ^nl , (3) 
n 

with /„ = Q-E„/kBT I ^-En/kBT normalized Boltzmann factor and the n-th eigen- 
state of the Hamiltonian. The eigenenergies {En\ are independent of time. The time depen- 
dence of the density matrix is through the vortex position dependence of the wavefunction. 
Then we have, to the lowest order in the vortex velocity, 

jV{t) = Tr l^/nV. ■ [IVro^n >< + >< Vro^n|]p| • (4) 

Here the operator Vro is the gradient with respect to the vortex position, and v„ = fq. Since 
the vortex position, or the position of the trapping potential, is the only reference point in 
the electron fluid system, we have the identity 

Vro*n({r,},ro) = - ^Vr,^'n({r,},ro) , (5) 

3 
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with Tj the position of the j-th electron in the system. Note that the total momentum 
operator is P = —ihj^j Vr^ = — «^Vr, Eq.(4) gives us 

jP{t) = ih 1^ /„[- < Vr^„|v„ • Vr^„ > + < ■ Vr^„|Vr^„ >]| 

= -thv, xEfnf n^'r,[VR X VR.M/„({r;.})M/:({r,})]|{r;}={r,} ■ (6) 

Following the same procedure as in Ref.p], we first reduce the N-body density to the one- 
body density matrix pi. 



~dt 

with 



P(t) = -ihw, X j d=^ri[Vr, X Vr;pi(r;;ri)]|r;=ri , (7) 



and the total electron number. Then we use the Stokes' theorem to evaluate the integral 



far away from the vortex core, and obtain 



jV{t) = Lh^j, X 7.jdY- ^[(Vr- Vr')pi(r';r)] 



L q^h^v^ X z , (9) 



which shows that the force felt by the electron fluid has the same magnitude as the Magnus 
force, but with the opposite sign. We will return to this point below. In the integration 
leading to Eq.(9) we have used the fact that the integrand is the momentum density. The 
two-fluid model has been employed to account for the fact that the momentum generated 
by the vortex corresponds to a supercurrent. The momentum flow calculation can also be 
performed near the vortex core by counting the flow of the energy spectrum, the so-called 
spectral flow method, and have been done in Ref. P|. The main idea is that, by linearizing the 
time-dependent Bogoliubov-de Gennes equation near the Fermi surface, a pseudorelativistic 
Dirac equation will be obtained. This linearization reduces the original 3-1-1 dimensional 
problem to an effective 1-1-1 dimensional one. As a vortex moves, there is a continuous flow 



of the energy spectrum, emerging from (or sinking into) the Fermi sea. The spectral flow rate 
is proportional to the vortex velocity. By counting their contributions to the momentum, 
the same result as Eq.(9) has been arrived at in Ref. |^. This procedure has been further 
verified by a model calculation.]^ The agreement between two seemingly totally different 
ways, looking from a safe distance away from the core and watching the spectral flow near the 
core, of calculating the momentum flow may first appear surprising. It has been, however, 
explicitly demonstrated in Ref. in a slightly different context in the ^He-A phase that they 
are indeed the complementary two ways of keeping track of the same physics, and they are 
equal as a result of the momentum conservation in the electron fluid system. In the present 
situation the link between those two ways has been discussed in Ref. 0]. It is the differential 
form of the momentum flow: 

d 

— (m*v^) + V/x = g„/iv^ X z 5^(r - ro) , (10) 

with the electrochemical potential as the sum of the chemical potential /xo, the fluid kinetic 
energy and the electric potential /i = /xo + ^^s+e*(0+^^ f dv' ■A{r' , t)), and the superfluid 
velocity distribution m*Vs = hVO — ^ A. Here m* and e* are effective mass and charge of a 
Cooper pair, respectively, and 6 is the phase of the superconducting condensate wavefunction. 
Eq.(lO) is gauge invariant, and is a precise statement about the phase slippage process 
due to a moving vortex. Its various consequences have been explored by Josephsonf^ and 
Anderson[^]. An identical equation for the special case of = in the neutral superfluid 
has been studied by Anderson]^]. We may call Eq.(lO) the Josephson- Anderson relation. 

It is clear now that a moving vortex feels a transverse force, the Magnus force; if one looks 
from the point of view of the electron fluid, the fluid also feels a force with the magnitude 
equal to but the sign opposite to the Magnus force. The two forces are acting on two different 
objects, the vortex and the electron fluid, and are the action-reaction forces. In the following 
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we strengthen this point by a straightforward proof. 



We write Eq.(2) in its equivalent form of the Ehrenfest theorem: 



d 



P(t) 



Tr[p{t)VnH] . 



(11) 



Using Vr-?/ = — Vro-f^, we obtain 



d 
di 



P(t) = TrlpmroH] 



(12) 



This shows that the two forces under discussion are indeed the action-reaction forces. In the 
last equality we have used the formal definition of the Magnus force0]. 

The present result, the equivalence between the Magnus force and the spectral flow as 
action-reaction forces, may seem obvious. Nevertheless it has never been explicitly spelled 
out in the literature. Instead, some recent work have treated the spectral flow as a way to 
cancel the Magnus force [§, ^ §, which is incorrect according to the present demonstration. 
It should be pointed out that the spectral flow is a counting of contributions from extended 
states. There is no involvement of the localized core states. This can be explicitly checked by 
expressing the one-body density matrix in Eq.(9) in terms of extended and localized states, 
and latter gives zero contribution, as having been noted in Ref. in the evaluation of the 
Magnus force. Incidentally, in Ref.|^ the Wess-Zumino term for a moving vortex has been 
identified as the same force as the one due to the spectral flow. Recent as well as earlier 
work have shown that the Wess-Zumino term gives exactly the Magnus force[|l^, 0, H, 
not the spectral flow. 

Comparing the Berry's phase calculation away from the vortex core with the spectral or 
momentum flow counting near the vortex core, we find that the former only depends on a 
few global properties of a superconductor, namely the topology of a vortex, and the latter 
is a rather detailed calculation. The topological constraints behave like conservations laws. 
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Results obtained under them should, and have to, be borne out by detailed calculations, 
which are concrete realizations. For the Magnus force, it is indeed the case. 
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